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Abstract 

AMS-02 provided the unprecedented statistics in the measurement of the positron fraction from 
cosmic rays. That may offer a unique opportunity to distinguish the positron spectrum coming 
from various dark matter (DM) annihilation channels, if DM is the source of this positron excess. 
Therefore, we consider the scenario that the DM can annihilate into leptonic, quark, and massive 
gauge boson channels simultaneously with floating branching ratios to test this hypothesis. We also 
study the impacts from MAX, MED, MIN, and DC diffusion models as well as from isothermal, 
NFW, and Einasto DM density profiles on our results. We found two parameter regions that can 
satisfy both AMS-02 e+ e + e - an d p/p datasets at 95% CL. i) Under the NFW-MIN combination 
with M x C [10,30] TeV. ii) Under the Einasto-DC combination with M x C [500,1500] GeV. 
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I. INTRODUCTION 


Up to now the evidences of the existence of dark matter (DM) only come from the 
gravitational effects. In the early 1950s, people noticed the anomaly in the rotational velocity 
in the Milky Way and other spiral galaxies, and one logical reason is that there are non- 
luminous matters suffusing our universe. Currently the Planck satellite [1] pinned down that 
26.8% of our universe is made up of DM. The Bullet cluster and N-body simulations may 
reveal some self-interaction of DM particles. 

From the particle physicists’ point of view, the DM particle might have interactions with 
the standard model (SM) particles, that will lead to interesting phenomenology. The DM 
direct detection experiments, like XENON100[2J, CDMS II[3] etc., can probe the interac¬ 
tions between the DM and nucleon. Searching for the DM signal as missing energies have 
been performed at collider experiments such as ATLAS ESE] and CMS [g]. In this work, we 
focus on DM indirect detection. Since the DM suffuses our galaxy, they probably annihilate 
with each other into SM particles. In this case, the antimatters, like positrons or antipro¬ 
tons, would be the ideal signals for searching for these processes, because the generation of 
antimatters is rare from ordinary astrophysical processes. The Alpha Magnetic Spectrom¬ 
eter (AMS-02) is a particle detector, which was designed to search for the antimatter and 
determine the composition and flux of the cosmic ray. The AMS-02 published the measure¬ 
ments of electron and positron results with unprecedented statistics in 2014. Their results 
confirmed the positron excess with energy higher than 10 GeV up to 500 GeV, which was 
first indicated by the PAMELA on 2008 [10], [ITj . On the contrary, based on PAMELA and 
most recent AMS-02 |9j observation there is no significant excess of antiproton in the cosmic 
rays from the prediction of conventional astrophysical processes. 

The DM annihilation is one of the viable interpretations for the observed positron excess, 
however it depends on which kinds of the SM particles that the DM particles annihilate into. 
The antiproton signal in general co-exists with the positron signal as parts of the annihilation 
products, that is severely constrained by the exisiting antiproton observations. This is one 
of the main difficulties in explaining the positron excess by DM annihilation. Furthermore, 
the newly released e ++ e - data from AMS-02 has so high statistics that it may be used to 
distinguish various DM annihilation channels. In this work, we set up a model-independent 
DM annihilation framework, in which the DM could annihilate into lepton channels (e + e _ , 
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p + ,p~, t + t~), quark channels (tt, bb), and massive gauge boson channels {W + W~, ZZ ) 
for a single annihilation with varying branching ratios (BR). Through fitting to the e ++ e - 
and p/p data, we investigate if this DM scenario can survive, if it does, which one is the 
preferred channel. Here we use the bb channel to represents other lighter quark channels, 
because with DM mass above few hundreds GeV, the positron and anitproton spectra from 
bb channel are similar to those spectra from other lighter quark channels. 

Note that AMS-02 very recently released the preliminary data on the ratio p/p [9], which 
also seems to have an unexpected rise at the first glance. However, it was shown to be 
consistent with the astrophysical background and conventional diffusion model within un¬ 
certainties [3QTf32] . We will include the most recent AMS-02 p/p data in our analysis. 

This work is organized as follows. We briefly describe the AMS-02 e+ e +e _ data and 
background models in the next section and DM model in Sec.Ill. The fitting results of e+ e +_ 
data shall be shown in Sec.IV. In Sec.V, we start considering the constraint from AMS-02 
p/p data. In Sec.VI, we discuss the influences from density profiles. Then we conclude in 
Sec.VII. 


II. AMS-02 POSITRON DATA AND BACKGROUND 

In September 2014, the AMS collaboration published their high statistics measurement 
of the g+ e _^ e _ in primary cosmic rays of 0.5-500 GeV [7J. The e ++ e - increases and exceeds 
the astrophysical cosmic ray expectation for energy larger than about 10 GeV. This trait is 
consistent with the earlier observations by PAMELA [10, Tl] and Fermi-LAT|12j. Further¬ 
more, the AMS-02 provides the most precise measurement so far such that it offers more 
information and tighter constraints on the origin of positron excess. 

The uncertainty of the high energy cosmic positron flux, which is originated from astro- 
physical sources and presumingly generated from the collision of high energy proton and 
interstellar medium (here we denote as the background), is still large. It is due to lack of 
knowledge of the sources and propagation models for the positron flux. The proton flux has 
been precisely measured by many experiments in the past and the most recently by AMS-02. 
The secondary electron and positron flux could be calculated straight-forwardly from the 
interactions of proton with interstellar material (mostly hydrogen and helium). However, 
there are large uncertainties in the propagation of the secondary particles from the origin of 
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production to the top of atmosphere (TOA). Difference diffusion models will give different 
positron spectra and as well as astrophysical background of e ++ e - when we deal with the 
AMS-02 positron fraction data. 

In this work, we will consider four different diffusion models, namely MIN, MED, 
MAX [22H23], and a diffusion convection model (DC) [25]. Each diffusion model is paired 
with a background, which is obtained by using the method mentioned above and computed 
through the Galprop package mm, and these backgrounds were shown in FigjTJ More 
description and the parameters for these models will be given in the next section. I 11 Fig. [lj 
we compare the four backgrounds of e+ e l from the four diffusion models. The ordering of 
background values of e ++ e - is MIN > MED > MAX, when the energy E e + > 10 GeV. For 
the DC case, the background is close to that of MED case but has a slight deficit, when 
E e+ > 100 GeV. 

The backgrounds of e +^ e - from MIN, MED, and MAX diffusion models are quite different 
from one another, in such a way that we anticipate that they are able to cover or represent 
larger portions of diffusion models in the market. Also, they are able to test whether our 
results are diffusion-model independent. The reason we adopt the DC diffusion model is that 
it utilized the more present B/C data to constrain the cosmic-ray propagation parameters. 


III. POSITRON SOURCE FROM DM ANNIHILATION 

The AMS-02 [7] provides a very precise measurement of the positron fraction, such that it 
could discriminate different DM annihilation channels and determine which one is preferred 
or ruled out. In order to do so, we consider the multi-annihilation channels DM scenario 
by assuming that the dark matter (DM) particles x could annihilate into seven standard 
model (SM) channels ti , bb, W + W~, ZZ , e + e _ , /i + /U, and t + t~ from a single annihilation 
process. The branching ratio for each channel can float between 0 and 1, as long as they 
satisfy the sum rule 

E br(//) = i. ( 1 ) 

f=t,b,W,Z,e,ii,T 

The relation between the total DM annihilation cross section and each channel is specified 
as 

( cn; >xx-U/ = BR (//) x ( av )tot ■ 
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These SM particles in the final state will further decay or hadronize into positrons, antipro¬ 
tons, or gamma-rays, which will be added on to the spectra obtained by the conventional 
high energy cosmic rays. Due to the distinctive spectral shape of these DM annihilation 
products, we may be able to separate these DM signals from the astrophysical cosmic rays 
and interpret as an indirect evidence of DM. In general, the leptonic channels e + e _ , p + p~, 
and t + t~ produce a relatively harder positron energy spectrum than the quark channels 
tt and bb. The positron spectra from massive vector boson channels W + W~ and ZZ sit 
between the leptonic and quark channels. 

We use the interface for model-independent analysis from micrOMEGAs [HI], where the 
Pytliia [TB] was embedded to generate the positron and antiproton spectra of each DM 
annihilation channel. Several DM density profiles will be considered in this work, including 
NFW [[26], isothermal and Einasto [ 281 f29| . Nevertheless, most of the time we focus on 
the NFW profile. The NFW and isothermal can be parametrized by 

( 0-7 )/a 


, \ [ r ol 7 ( 1 + ( r &/ r s) a ) 

PM = Po{-} { ! + (r/r j„ } 


where p Q = 0.3 [GeV/cm 3 ] is the DM density at the Sun, r 0 = 8.5 [kpc] is the distance from 
the Sun to the Galactic Center (GC). The NFW (isothermal) density profile can be obtained 
by choosing the parameters a = l,/3 = 3,7 = l, and r s = 20[kpc] (a — 2, /3 — 2,j — 0, and 
r s = 3.5[kpc]). On the other hand, the expression for the Einasto density profile is 


P(r) 


P . e {-mr-A} 


with p s = 0.3 [GeV/cm 3 ], r s = 8.5 [kpc], and a = 0.17, that will give p(r = r 0 ) = p & . 

In micrOMEGAs [15], the energy spectrum of the positrons is obtained by solving the 
diffusion-loss equation while keeping only the two dominant contributions: spatial diffusion 
and energy losses, 

- V • (K(E)VM - (b(E)M = Q e+ (x,F), (2) 


and 

K(E) = K 0 (3(E)(R/lGeV ) s , 


where (3 and R = p/q are the velocity and rigidity of the particle, respectively. The positron 
energy loss rate b(E) is dominated by the synchrotron radiation in the galactic magnetic 
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field and by inverse Compton scattering on stellar and CMB photons. The diffusion zone 
of cosmic rays is represented by a cylinder with thickness 2 L and radius R (here we use 
R — 20 [kpc]) [22]. The term Q e + on the right hand side of Eq. [2] is the positron source term. 

The three diffusion models (MAX, MED, and MIN) (221TM] that we are going to consider 
in this work, correspond to the diffusion parameters: 

MAX : L = 15 , [kpc], K 0 = 0.0765 , [kpc 2 /Myr], 5 = 0.46 , 

MED : L = 4 , [kpc], K 0 = 0.0112 , [kpc 2 /Myr], 6 = 0.70 , 

MIN : L = 1, [kpc] , K 0 = 0.0016 , [kpc 2 /Myr] , 5 = 0.85 , 


all of them satisfy the B/C observations. Since the AMS-02 also provided the measurement 
of the B/C ratio, for more comprehensive discussion of the influence of diffusion models, we 
add one more diffusion model, the DC model from Ref. [25], where the diffusion parameters 
were determined by fitting to the most recent AMS-02 B/C data. 

Under the DM with the multi-annihilation-channel scenario, it will contribute to the 
source term Q e + in Eq. (J2]) through the expression 


Qe + M = V 


Pdm 

mT 


Mm E BR (//) 

f=t,b,W,Z,e,ii,T (l£J e + 


(3) 


where dN^ + /dE e + is the positron energy spectrum from each annihilation event XX //, 
and rj = 1/2 (1/4) for (non-)identical initial state. The p^M is the DM density profile, and 
M x is the DM mass. 


IV. AMS-02 POSITRON FRACTION DATA AND FITTING RESULT 

Here we attempt to use the DM scenario with multi-annihilation channels, each of which 
has a varying branching ratio subject to the sum rule in Eq. (Jlj), to fit the AMS-02 e ++ e - 
data [7j. The varying parameters are the DM mass, the annihilation cross section, and the 
branching ratios of seven annihilation channels: M x , (crv) tot , and BR(//). However, the 
branching ratios should satisfy the sum rule in Eq. (Jl|, therefore, overall we have 8 free 
parameters. 

First, we fix the DM density to the NFW profile and use various diffusion models (MIN, 
MED, MAX, and DC) to compute the positron flux at TOA from DM annihilation. The 
positron source term of DM annihilation channels can be described by Eq. 0 . Then we 
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feed the source term into the diffusion equation (|2]) , and it is solved numerically by the 
rnicrOMEGAs. Finally, we stack the DM-produced electron and positron flux on top of 
each diffusion-corresponding background flux to obtain the values of e ++ e - ■ Then we are 
able to compare with the AMS-02 positron fraction data. 

We use the \ 2 and p-value as our statistical measures to describe the deviation of our 
model from the data and the goodness of the fit, respectively. [] Our scanning strategy 
is fixing certain values of M x and (crv) tot , meanwhile allowing the branching ratio of each 
channel to float as long as they satisfy the sum rule to obtain the minimal y 2 or maximal 
p-value. Repeating the whole procedures to scan the parameter space of M x and (crv) tot . 
The results are shown in Figs. mm and [5j corresponding to MED, MAX, MIN, and DC 
diffusion models, respectively. 

In the panel of (M x , ( av )) in each of the Figs. m and [5j we show the 68.3% (blue, 
0.32 < p — value < 1.0) 95% (green, 0.05 < p — value < 0.317), and 99% (yellow, 0.01 < 
p — value < 0.05) as regions of confidence level (CL) of the best-fitting points with varying 
branching ratios of the DM annihilation channels. In these panels, independent of the 
diffusion models, they show the narrow window in (crv) tot restricted by AMS-02 e ++ e - data 
for each value of M x , despite of the variation of seven annihilation channels. Also, the total 
annihilation cross section (crv) tot strongly correlates with and is roughly proportional to the 
M x , but it is much larger than the thermal DM cross section, which is about (w) thermal — 
1 pb ~ 3 x 10~ 26 cm 3 /s. 

In the panel of (M x , Br) in each of the Figs. 2pp , and[5j the light blue-shaded area shows 
the p-values of the best-fit point for the corresponding M x , and the p— value is labelled on 
the right-hand y-axis of the panel. It manifests itself in the (M x , Br) panel that once the 
DM mass is below 300 GeV, the p- value is always less than 0.05, independent of the diffusion 
models, and no matter how the combination of the annihilation channels changes. In other 
words, we can rule out the DM mass below 300 GeV as an explanation for the AMS-02 
e ++ e - excess at 95% CL. This is because the last two AMS-02 e+ e _^ e _ data points with the 
highest energy between 260 GeV to 500 GeV cannot be explained by DM annihilation for 
1 Assuming the goodness-of-fit statistics follows y 2 statistics, the p-value for the hypothesis is given by m 

/»oo 

P= f(z-,n)dz 


where n is the degrees of freedom and 


f(z;n) = 


z n/2-l e ~z/2 

2 n / 2 r(n/2) ' 
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the DM mass below 300 GeV, due to the cliff shape of the positron energy spectrum from 
DM annihilation around M x . On the other hand, the upper limit for the DM mass, which 
are diffusion model dependent, are 5.5, 60, 50, and 13 TeV for MAX, MED, MIN, and DC, 
respectively. 

There are two panels of (M x , Br) in each of the Figs. m3 and [5j one of them shows 
the branching ratio of individual channel for the best-fit point at the corresponding M x , 
while the other one shows the branching-ratio sum over the leptons(e + e _ + p + /i _ + t + t _ ), 
quarks (bb+tt), or massive gau.ge-boson{W + W~+ZZ) channels. When the DM mass is below 
10 TeV, the leptonic channels dominate the branching ratio (BR(leptons) ^ 60%), except 
for the region of M x < 1 TeV in the MAX diffusion model, where the quark channels are 
able to compete with leptonic channels. Nonetheless, among these four diffusion models, the 
quark and gauge-boson channels still play non-negligible roles because (BR(gauge bosons) + 
BR(quarks) > 20%). Since (&v) tot ~ O (10 3 ~ 10 4 pb), even 10% of the W + W~, ZZ, 
and bb channels would generate enormous antiproton flux, and that is why it is severely 
constrained by the antiproton data na. Therefore, the largest p-value point for each M x 
would be ruled out by the antiproton data. Among these four diffusion models, only the MIN 
diffusion model at the region with M x < 350 GeV, where the three leptonical channels can 


fit the AMS-02 


e++e~ 


with p —value larger than 0.05. Another feature is that once the DM 


mass goes above 30 TeV, the tt channel becomes more and more important and eventually 
becomes the major channel, because of the energy spectrum of tt is relatively softer than 
the other channels. 

In the following, we consider how the AMS-02 e+ %_ data can distinguish among different 
channels. The panel of (M x , y 2 ) in each of the Figs. mi and [ 5 ] shows the chi-squares 
minimum with respect to M x if the DM can only annihilate into a single channel. The 
solid lines represent the single-annihilation channel scenario, while the dashed line shows 
the chi-square minimum of multi-annihilation channel scenario. We could see that the 
DM with multi-annihilation channel scenario can reduce the y 2 significantly, although more 
parameters are involved. When the DM mass is less than 1 TeV, the t + t~ single-annihilation 
channel scenario can provide a relatively smaller value of y 2 than the other channels, and 
the minimum y 2 occurs at around M x =650, 800, 500, and 450 GeV for MAX, MED, MIN, 
and DC diffusion models, respectively. The two massive vector-boson channels attain their 
y 2 minimum at about M x = 7, 10, 10, and 5 TeV for MAX, MED, MIN, and DC diffusion 






models, respectively, because they generate moderate positron energy spectra in between the 
leptonic and quark channels. The ti channel, on the other hand, generates a further softer 
positron energy spectrum, such that the y 2 minimum occurs at around M x = 25, 40, 35, and 
20 TeV, for the MAX, MED, MIN, and DC diffusion models, respectively. Once M x > 20 
TeV, the curve of ti is in unison with the blue-dash curve, and the ti channel takes the full 
responsibility for fitting the data, and other channels become irrelevant. From these results, 
we can see that the lepton channels are well separated from quark and massive gauge boson 
channels, and therefore we anticipate that the AMS-02 e ++ e - data is able to distinguish the 
leptonic channels from the other channels in the DM multi-annihilation channel scenario. 
However, the y 2 minima of W + W~ and ZZ are almost overlapping, hence it is hard to 
distinguish them from each other. In the bottom-left panel of Fig. a with 300 GeV M x % 
450 GeV, pure leptonic channels are solutions of AMS-02 positron fraction data under MIN 
diffusion model and NFW profile. 

For MIN diffusion model in Fig. |4j there is no contradiction between the panels of (M x , 
Br) and (M x , y 2 ) around M x = 12 TeV. At the first glance, the panel of (M x , y 2 ) points 
out that the WW and ZZ channels are the most important ones to fit the data at this M x . 
However, the Br(e + e“) is the largest from the panel of (M x , Br). The reason is that the 
e + e~ channel almost decouples at this DM mass region, and so the positron spectrum from 
it is way above the energy scale of AMS-02 e ++ e ~ data, such that dramatically change of 
(av) x ^ e+e - does not influence the fitting too much. Therefore, the WW and ZZ channels 
are still the major channels to fit the data around M x = 12 TeV. 

The panel of (Br (l + l~), p) in Fig. [2]shows how many percentages of the leptonic channels 
are required to fit to the AMS-02 e ++ e - data with the corresponding M x in the DM multi¬ 
annihilation scenario, while marginalizing all the massive gauge boson (VV) and quark (qq) 
channels. Take M x = 1000 GeV, for example, the y 2 minimum occurs with Br(/ + Z~) ps 
75±23% and Br(qq) + Br(VV) « 25%^% at 95% CL in the multi-annihilation scenario, while 
the single-annihilation channel scenario is the r + r^ channel but at a much larger y 2 . Two 
inferences we can draw from here for M x — 1000 GeV: (i) the major channels are the leptonic, 
and (ii) the three lepton channels alone, Br(/ + / - ) = 100%, could not fit well to the data. 
Instead, they need help from softer positron spectra of quarks or gauge bosons. When M x < 
10 TeV, the three leptonic channels dominate. On the other hand, for M x > 10 TeV, the 
VV and qq channels gradually become the major contributions. For instance, when M x = 
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30 TeV, Br(Z + Z - ) « 70^]o%. Yet, without any lepton channels the VV and qq still have 
p-value = 0.3. In another words, the AMS-02 e ++ e - data constrain more severely on the 
branching ratio of leptonic channels for M x < 10 TeV. 

In Fig. 6 we would like to demonstrate how different background of p+ e + e _ originating from 
different diffusion models will modify the conclusions that we made in previous paragraph. 
For all three diffusion models MAX, MED, and MIN at M x ~ 1 TeV, the values of Br(/ + /“) 
are still confined in a narrower window with about ±20% deviation at 95% CL from the 
best-fit value of Br(7"% _ ). On the other hand, for M x > 10 TeV the confinement is much 
loose such that the values of Br(7"% - ) can go down to zero in MED and MIN cases. Another 
feature we can see from this figure is that in the MIN (MAX) diffusion model, the values of 
Br(/ + Z - ) concentrate at slightly larger (smaller) values than those of the MED model. 


V. AMS-02 ANTIPROTON-PROTON RATIO DATA AND FITTING RESULTS 

During 15-17 April 2015 at CERN Geneva, the AMS collaboration presented the pre¬ 
liminary results on the spectrum of the antiproton-to-proton ratio for proton kinetic energy 
(KE) up to 450 GeV 0, with much higher precision than other measurements in the past. 
The data points with KE between 20 and 450 GeV show a flat structure. Comparing the 
AMS-02 and PAMELA p/p data, with KE larger than 5 GeV up to 180 GeV, these two 
datasets are consistent within the errors. 

The derivation of the antiproton background is analogous to the positron case. First, 
we specify the diffusion model and determine the proton flux from the data, then compute 
the antiproton flux based on Galprop [13] Q3]. Various backgrounds originating from the 
corresponding diffusion models are shown in Fig. [7j There is a moderate excess of the AMS- 
02 p/p data above various backgrounds. Comparing among these four backgrounds, the one 
from the DC diffusion model provides a slight better fit to the AMS-02 p/p data. 

We keep four annihilation channels (ft, bb, W + W~, and ZZ) of DM for the the calcu¬ 
lation of antiproton ratio, because these channels can generate antiprotons after decay and 
hadronization, while the e + e _ and p + p~ channels are irrelevant and the contribution from 
the r+r' channel is negligible. We therefore define (. crv) p/p = 'Eff=ti,bb,ww,zz(w) xx ^ff- 
We take this multi-annihilation DM scenario to fit to the AMS-02 p/p data, and use only 
the data with KE larger than 5 GeV in the fitting. The results in terms of CL regions in the 
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(M x , (crv)) panels with MAX, MED, MIN, and DC diffusion models are shown in Fig. [8j 

In Fig. |8j we can see a strong and almost linear correlation between M x and (c r v) p ^ p 
in all four panels. The 95% CL lower and upper limits on M x depends on the diffusion- 
model. The 95% lower limit on M x = 850, 600, 1000, and 450 GeV for MAX, MED, MIN, 
and DC diffusion models, respectively. The value of (av) P//p required to fit the AMS-02 
p/p data is always smaller than the ( crv) tot required to fit the AMS-02 e+ e + e _ data with the 
corresponding M x , except for the MIN diffusion model with M x > 10 TeV. Nevertheless, the 
values of ( crv) tot and (av) p / p are of the same order of magnitude. The values of (crv) p , p for 
the MIN diffusion model are about two orders of magnitude larger than the other models, 
mainly because of the fact that the diffusion zone of MIN is thinner than the other diffusion 
models such that particles are easier to escape the diffusion zone during propagation. It then 
requires larger values of (crv) p , p to compensate for the lost particles. Thus, for M x > 10 
TeV the value of (<jv) pPp for the MIN diffusion model is of the same order as the value of 
(av) tot in Fig. [4j In other words, in the NFW-MIN case for M x > 10 TeV, it can tolerate 
large branching ratios of massive vector boson and quark channels without conflict with the 
upper bound from AMS-02 p/p data. 

The 95% CL upper limit on (av) p j p shown in Fig. |8]are also diffusion-model dependent. 
Taking the MED diffusion model as an example with M x = 1(10) TeV, the upper limit on 
( <Jv) p / p ~ 20 (1000) pb, which is about 1/30 of the favored value of (crv) tot in the AMS- 
02 e+ e + c - case. In other words, if we try to fit both AMS-02 p/p and e+ e + e _ data, the 
BR (qq) + BR(W) < 3% or Br(Z + /“) > 97% is required. However the panel of (Br(Z + Z — ), p) 
in Fig. [2] tells us that the parameter region satisfying the above requirement is very narrow: 
only when M x < 1 TeV and the p —value is less than 0.1. After scrutinizing the parameter 
space of the multi-annihilation channel DM scenario, it cannot provide simultaneously good 
fits to both the AMS-02 e + +e - and p/p data with a p —value larger than 0.05 based on the 


NFW-MED combination. 

A simultaneous solution for both AMS-02 e +l e - and p/p data could be found under the 
NFW-MIN combination for M x C [10,30] TeV, as shown by the compensating color in the 
panel of (M x , (crv)) in Fig. [dj The upper two panels in Fig. 11 show the a benchmark point 
of the solution for both AMS-02 e+ e + e - and p/p data, with both p-values larger than 0.05: 
M x = 23 TeV with composition of 22 %t + t“, 39% tt, and 39% ZZ channels. 

After the comprehensive parameter survey of the multi-annihilation channel DM, we did 
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not find any solution to simultaneously explain the AMS-02 e+ e + e _ and p/p data under the 
case of the NFW profile combined with MAX, MED, or DC diffusion models. 


VI. DIFFERENT DARK MATTER DENSITY PROFILES 


So far in the above sections, we only focus on the NFW DM density profile. In this 
section, we investigate the consequences of using different density profiles by adopting the 
isothermal and Einasto profiles. 

Among these prohles, we have normalized the DM density to p(r@ = 8.5 kpc) = 
0.3 GeVcm -3 at our solar system, where r Q is the distance of our solar system from the GC. 
Note that the largest uncertainty among various prohles happens within 1 kpc around the 
GC. Around the vicinity of our solar system, the isothermal and NFW prohles are almost 
identical, but however some deviations between Einasto and NFW happens. Hence, the 
positron spectra from DM annihilation from isothermal and NFW prohles are quite similar, 
but different from the Einasto prohle. In general, the Einasto prohle will give a softer 
positron spectrum than the NFW and isothermal. This feature could make it easier for the 


multi-annihilation channel DM scenario to ht simultaneously to both AMS-02 


e++e~ 


and 


p/p data, because the missing ingredient in the NFW prohle is that the positron spectrum 
from leptonic channels is not soft enough for the e ++ e - data. In the following, we take the 
combination of the Einasto prohle and the DC diffusion model to demonstrate this feature, 
since the DC model utilizes a more recent B/C observation to constrain their cosmic-ray 
propagation parameters. 

Based on the Einasto-DC combination, the AMS-02 e+ e +_ and p/p fitting results are 
shown in Figs. [9] and [lOj respectively. The benchmark point of multi-annihilation channel 


DM that can explain both data sets is illustrated in the lower two panels of Fig. [TTJ Com¬ 
paring the Fig. [5] (NFW) with Fig. [ 9 ] (Einasto), larger values of M x > 10 TeV are allowed 
at 95% CL in the Einasto prohle. From the panels of (M x , y 2 ), the fitting of the single 
channels of t + t~, VV and qq are all improved in the Einasto prohle than those in NFW 
one. Also, because of the improvement in VV and qq channels, the region of M x > 10 TeV 
is allowed. Another difference in the panels of (M x , y 2 ) between these two prohles is that 
the corresponding values of M x of the minimum of each channel shifts to larger values in 
the Einasto case. This comes from the fact that the Einasto prohle makes positron spectra 
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from each channels softer. In the panels of (M x , Br) in Fig. [9j when M x less than 1200 
GeV, the Br(Z + Z - ) of the best-fit point can be larger than 95%, and so Br(W) + Br (qq) is 
suppressed below 5%. This region provides a nice habitat to accommodate the AMS-02 p/p 
data. In the panel of (M x , (crv)), the red long-elliptical contour indicates the allowed DM 
mass and annihilation cross section that can explain well these two AMS-02 observations. 
The DM mass range is M x C [500,1500] GeV. 


In Fig. 10, we show the results of fitting to the AMS-02 p/p data under the Einasto-DC 
combination. Comparing the panel of (M x , (crv)) in this figure to the bottom-right panel 
(corresponding to the NFW-DC combination) in Fig. [8j the 95% CL allowed range for DM 
mass does not change much, 400 GeV ^ M x < 15 TeV. However, because of the DM density 
p(r) of the Einasto profile is larger than that of NFW for r < r Q , the entire CL region 
shifts down by a factor of 2 for the Einasto case, indicating that the value of ( <Jv) tot for the 
Einasto case is half of that for NFW. Similar observation can be seen from the fitting to 
the AMS-02 e ++ e - data, by comparing the panels of (M x , (crv)) in Figs. [H] and 9 Another 
observation from the two panels of (M x , Br) in Fig. 10 for the best-fit point is that when 
M x < 1700 GeV, the massive gauge boson channels take the most responsibility to fit the 
data. However, when M x > 1700 GeV, the quark channels take over the major role to do 
the fitting. 

In summary, the two major differences between the Einasto and NFW profiles are (i) the 
former makes the positron spectra from DM annihilation softer than the latter, and (ii) in 
order to produce the same amount of positron or antiproton flux from DM annihilation, a 
slightly smaller annihilation cross section is required in the Einasto profile. 


VII. CONCLUSIONS 

In this work, we have postulated a DM scenario with multi-annihilation channels to 
account for the positron fraction excess and antiproton-to-proton ratio observed in cosmic 
rays. By using the most updated e ++ e - and p/p data from AMS-02, we have performed 
a comprehensive analysis to determine the preferred DM annihilation channel in terms of 
statistic quantities y 2 and p-values. We have considered the DM annihilation into tt, bb, 
W + W~, ZZ, e + e~, p + , and t + t~ channels from a single annihilation process, and taken 
the branching ratio of each channel and the total cross section as the fitting parameters. We 
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focus on two aspects: i) Which channel is preferred by the AMS-02 e +l e - data, and how 
strong the confinement of branching ratios can be implied, ii) Under this multi-annihilation 
channel DM scenario, we are looking for a simultaneous solution for both AMS-02 e+ e _^ e _ data 
and p/p data. Several diffusion models (MAX, MED, MIN, and DC) and DM density profiles 
(isothermal, NFW, Einasto) are invstigated in our analysis. Uncertainties in astrophysical 
backgrounds are obtained via combinations of various density profiles and diffusion models. 

Our findings are summarized as follows. 

1. The allowed region of DM mass M x and annihilation cross sections (av) tot that can 
provide good fits with p-value > 0.05 to the AMS-02 e+ e + e _ data is very broad. The 
range of M x extends from 300 GeV to a few tens of TeV, and {crv) tot extends from 
0(100) pb to O(10 5 ) pb. However, there is strong and almost linear correlation between 
the M x and (crv) tot . 

2. The DM mass below 300 GeV is ruled out at 95% CL by the e ++ e - data, and this lower 
limit of M x is diffusion-model and density-profile independent. On the other hand, 
the upper limit for the DM mass strongly depends on the diffusion models and density 
profiles. Pure leptonic channels are solutions of AMS-02 positron fraction data under 
combination of NFW-MIN with M x C [300, 450] GeV or combination of Einasto-DC 
with M x C [300, 1800] GeV. 

3. The high statistical measurement from AMS-02 has the capacity to distinguish among 
different annihilation channels. For a few hundred GeV to several TeV DM, the r + r“ 
channel is the ideal one to fit to the e+ e +_ data. On the other hand, for several TeV 
to 10 TeV DM, the W + W~ and ZZ are the favorites. Beyond 10 TeV the bb and tt 
channels begin to take over and completely dominate very quickly. 

4. Based on the multi-annihilation channel DM scenario, for all MAX, MED, and MIN 
diffusion models, when M v ~ 1 TeV, the AMS-02 _ data has the best confinement 
of the branching ratio of the sum of three leptonic channels, and Br(/ + /“) is within 
±20% deviation from the best-fit value at 95% CL. However, the value of Br(/ + V) for 
the best-fit point depends on the diffusion model. 

5. The DM scenario with multi annihilation channels under the combinations of the 
MAX, MED and DC diffusion models, and the NFW density profile cannot satisfy 
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both AMS-02 e + e + e - and p/p data at 95% CL. However, under the NFW-MIN and 
Einasto-DC combinations, we find windows in, respectively, M x C [10,30] TeV and 
M x C [500,1500] GeV that can fit well with p-value > 0.05 to both datasets, For the 
Einasto-DC case, the t + t~ channel dominates and the massive gauge boson and quark 
channels are suppressed. On the other hand, for the NFW-MIN case the ZZ and tt 
are the two major channels. 

During the last stage of the manuscript, we came across a few similar works in model- 
independent approach [ 3011331435] and in specific models [36, 37]. 
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FIG. 1. The astrophysical backgrounds for the observable of <3? e +/(<i> e - + <P e +) using different 
diffusion models, corresponding to the MIN, MED, MAX, and DC diffusion models. Red crosses 
are the AMS-02 e _ e + e+ data m and blue crosses are the PAMELA e —_|_ e + data naan. 
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FIG. 2. Fitting results of the multi-annihilation channel DM scenario to AMS-02 data 

for the combination of NFW-MED. Description goes from top to bottom and from left to right, 
i) Upper-left: the allowed region in the plane of (M x , (av))tot- The blue region indicates 0.32 < 
p —value < 1.0 (68.3% CL), the green region 0.05 < p —value < 0.317 (95% CL), and the yellow 
region 0.01 < p — value < 0.05 (99% CL), ii) Upper-right: (M x , BR). The branching ratio of each 
channel for the best-fit points, labelled on the left y-axis. The blue-shaded region is the p —value 
of the best-fit points labelled on the right y-axis. iii) Middle-left: (M x , Br). Same as ii), but for 
summing over the leptonic, quark, and massive gauge boson channels, iv) Middle-right: (M x , % 2 ). 
Each of the solid curves shows the x 2 i R the fitting to the AMS-02 e -+ e + data using only a single 
DM annihilation channel, namely, bb. tt, e~e + , p~p + , t~t + , ZZ, or W~W + . The dashed curve 
represents the x 2 from the multi-annihilation channel DM scenario, v) Lower: (Br(Z _ Z + ), p). The 
sum of leptonic branching ratios Br(Z“Z + ) = Zl/= e u,r Br(//) vs the p-values for a number of M x 
values, labelled right next to the curves while marginalizing all other fitting parameters. 
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FIG. 3. The same as Fig. [2] but for the combination of of NFW-MAX, except without the last 
panel. 
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FIG. 4. The same as Fig. [2] but for the combination of of NFW-MIN, except without the last 
panel. The compensating color in upper-left panel in the plane of (M x , (av }) shows the allowed 
region of fitting simultaneously to both AMS-02 e ~+ e + and p/p data with p-value > 0.05 for each 
dataset. 






FIG. 5. The same as Fig.[2]but for the combination of of NFW-DC, except without the last panel. 
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FIG. 6. The same description as the lower panel in Fig. [2j except for using all different diffusion 
models. MIN: dotted curves; MED: solid curves; MAX: dashed curves. The color bar on the 


right-hand side indicates the DM mass M x . 



FIG. 7. Astrophysical backgrounds for the observable p/p data with the MIN, MED, MAX, and 
DC diffusion models. The red dots are the AMS-02 p/p preliminary data. The blue dots are the 


PAMELA p/p data. 
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FIG. 8. The fitting results in the plane of (M x , {crv)) of the multi-annihilation channel DM 
scenario to the AMS-02 p/p data. The relevant channels are tt, bb, W + W ~, and ZZ\ (av) = 
Ylf=t,b,w,z ( av ) X x-^ff' The blue re §i° n indicates 0.32 < p —value < 1.0 (68.3% CL), the green 
region 0.05 < p —value < 0.317 (95% CL), and the yellow region 0.01 < p —value < 0.05 (99% 
CL). The panels are: upper-left for the combination of NFW-MAX, upper-right for NFW-MED, 
bottom-left for NFW-MIN, and bottom-right for NFW-DC. 
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FIG. 9. The same as Fig. [2] but for the combination of of Einasto-DC, except without the last 
panel. The red region in the upper-left panel in the plane of (M x , (av)) shows the allowed region 
of simultaneously fitting to both AMS-02 e _ e + e+ and p/p datasets with p-value > 0.05 for each 
dataset. A zoom is displaced in the upper-left corner of that panel. 
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FIG. 10. Fitting results using the combination of Einasto-DC with respect to the AMS-02 p/p 
data. Upper-left: (M x , (av )) - confidence-level regions, labeled in the same way as Fig. [8j Upper- 
right: (M x , Br) - fitting to the AMS-02 p/p data alone with tt, bb, W + W~, and ZZ cha nn els. 
Lower: (M x , Br) - similar to the upper-right one, but summing over the similar channels: quark 
= E and massive gauge boson = J2f=w,z Br (//)- 
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FIG. 11. The left-handed panels are for the fitting results to the AMS-02: data, while the 

right-handed panels are for the fitting results to the AMS-02 p/p data. Upper two panels: under the 
NFW-MIN combination with a benchmark point within the compensating color region in the ( M x , 


(av)) panel in Fig. 


Br (ft) = 39% 


M x = 23 TeY, (av) tot = 1.330 x 10 5 pb, Br (t+t") = 22%, Br (ZZ) = 39%, 


Lower two panels: under the Einasto-MED combination with a benchmark point 


within the red contour in the (M x , {av )) panel in Fig. 


M x = 800 GeV, (av) tot = 5.962 x 10 2 


pb, Br(e+e-) = 10%, Br(p+p~) = 10%, Br(r+r-) = 78%, Bt(W + W~) = 2% 
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